Percolation Effects in Very High Energy Cosmic Rays 
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Most QCD models of high energy colhsions predict that the inelasticity K is an increasing function 
of the energy. We argue that, due to percolation of strings, this behaviour will change and, at 
^/s ~ 10* GeV, the inelasticity will start to decrease with the energy. This has straightforward 
consequences in high energy cosmic ray physics; 1) the relative depth of the shower maximum X 
grows faster with energy above the knee; 2) the energy measurements of ground array experiments 
at GZK energies could be overestimated. 
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Most QCD-inspired models of multiparticlc production 
predict, in hadron-hadron and nucleus- nucleus collisions 
at high energy, an increase with energy of the inelastic- 
ity parameter, K = I — xp, where xp is the momen- 
tum fraction carried by the fastest particle: Multiple 
Scattering Model Dual and string Models and 
Minijet Model Recently, several papers appeared 
studying collective and non-linear QCD effects, based on 
the Colour Glass Condensate Model (jj, Reggeon Cal- 
culus 5|, and Strong Field string Model 0, These 
models predict large stopping power and a decrease of 
the momentum fraction carried by fast particles. 

Two aspects are essential in the String Percolation 
Model 11,01 that we use here: 1) At low energy (or den- 
sity) leading valence strings are produced. As the energy 
increases, energy is drawn from the valence strings to 
produce, centrally in rapidity, sea strings. As in all the 
models mentioned above, the inelasticity increases with 
the energy. 2) At very high energy, above the percolation 
threshold, percolation leads to the formation of a large 
cluster of strings and to the production of faster parti- 
cles. As a consequence the inelasticity starts to decrease 
with the energy. 

While the relatively low energy regime, with K in- 
creasing, is similar to the models already mentioned, the 
higher energy regime, with decreasing K and the regen- 
eration of the fast particles, is new and has some straight- 
forward consequences in cosmic ray physics. 

In the String Percolation Model ^ for hadron- 
hadron collisions, at low energy, valence strings are 
formed, forward and backward in the centre-of-mass, 
along the collision axis, containing most of the collision 
energy. As the energy increases, additional sea strings, 
central in rapidity, are created, taking away part of the 
energy carried the valence strings. In the impact param- 
eter plane all the strings look like disks, and we have to 
deal with a two dimension percolation problem. 

The relevant parameter in percolation theory is the 



transverse density, rj [l£ 
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where r is the transverse radius of the string, R the ef- 
fective radius of the interaction area. TV^, the average 
number of strings, depends on the density (centrality) 
and on the energy. The strings may overlap in the inter- 
action area, forming clusters of N strings. If ?y ^ 1, the 
average number of strings per cluster is < >~ 1. If 
77 > 1,< TV >~ iVs. 

If n is the particle density for one string, TnT the av- 
erage transverse mass produced from a single string, and 
there are N s strings, one expects: 
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with a colour summation reduction factor |llL Il2| , 
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The particle density does not increase as fast as A''^ (this 
corresponds to the saturation phenomenon 0), and < 
rriT > slowly increases with energy and density. These 
features are seen in data (see, for instance, (isf). 

Following 01 , let us consider proton-proton collisions 
and write for the invariant s, 



5= (P1+P2)' 
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where Pi, 2 are the momenta of the protons, P =| Pi | = | 
P2 \, m is the proton mass, and Ay the length of the 
rapidity "plateau". For a string made up of two par- 
tons with Feynman-a; values X- and x^, we have, for the 
centre-of-mass energy of the two partons, si — x^x^s. 
Assuming for simplicity a symmetrical situation around 
the centre-of-mass, x- ~ x+ — x, the string centre-of- 
mass energy is si = x'^s, and we can write the length of 
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the rapidity plateau for the string as: 
Ayi = Ay + 2 lux. 



(5) 



If we write x = 2p/^, where p is the absolute value of 
the momentum of the partons, we obtain: 
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We shall assume, for sea strings, that p — constant, which 
implies a; ~ 1/ ^/s, in agreement with the rise of parton 
density at small x. Note that for valence strings, before 
energy degradation, p ^ P, and the full phase space is 
occupied by the valence strings. 

If strings overlap in the interaction region, and if < 
TV > is the average number of strings per cluster we have, 
generalising 



Ay 



<N> 



Ayi + 2 In < > 



(7) 



At low energy/density < >~ 1 and only short strings 
are formed, not contributing to cosmic ray cascades. At 
high energy/density < N >~ Ng, percolation occurs and 
the situation changes. 

In percolation theory the average number < iV > of 
strings per cluster is related to the average area < A >, 
in units of r^, occupied by a cluster |f^. 



< N >=< A > 
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and rjc — 1.15 is the transition point and a ~ 0,85 is a 
parameter controlling the slope of the curve at the tran- 
sition point, with /(ry) changing from to 1 at ~ Ty^- 
We note that when rj 0,< A >~ 1, and when 
77 — + oo, < v4 >~ {R/rY . This kind of parameterisation 
was tested in ;15j . 

The energy, in the centre-of-mass, carried by the pro- 
duced particles from sea strings is given by: 



with < A> given by 0| 

<A>= /(t?) 
where /(jy) is a percolation function, 
f{ri) = (1 + e-("-''<=)/'^ 
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and we obtain, making use of ^ and subtracting the 2 
valence strings. 
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If we now require that asymptotically all the energy is 
carried by the percolating strings. 
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we obtain, from (|3J|, Q and ((T^ . 
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with A = 2/7. 
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As Ns is proportional to the high energy bare Pomeron, 
the value of the intercept ap is related to A: — 1 = 
A. This result is consistent with results from the Colour 
Glass Condensate Model ,17j. 

In order to have an estimate of the inelasticity K we 
make use of the idea that, in the fragmentation of the 
string, produced particles are ordered in decreasing ra- 
pidity, and the fraction of momentum carried, relative to 
the momentum left, is always the same [l^ . 

At small v^, when the valence strings carry all the 
energy, the fastest particle {F) is the leading particle (L) 
and xp = XL = a = const., with < a < 1. When sea 
strings are produced, carrying an energy Ecai, we have. 
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When the strings percolate, Ecm ^ 2 | P | and for the 
fastest percolating particle (P) we have: 
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As Ecm is an increasing function of y^, xl decreases 
with the energy and xp increases with energy. Thus, 
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1 — XL; for XL > Xp 
1 — .Tp , for XL < Xp 
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In order to implement the model (equations |(2J), Q, 
(|12|l and IjlSI) to H17|l ). we have to establish the relation 
between Ng and 77(eq. iQ)), and to fix the parameters 
r/i?, rriT and Ayi. At some low energy threshold, ~ 
10 GeV, we have just the valence strings and N s = 2. At 
00, Ns ~ with A given by equation 114|l . We 
then write: 
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St 



(18) 



where the parameter 6 = 1.37 was adjusted to agree with 
the data on dn/dy. The remaining parameters were fixed 
to reasonable values: r/R = 0.2, rriT = 0.78 and Ayi — 
6. In this way, H13|l was exactl^atisfied. 

In Fig. n] the dn/dy data 19] is compared with the 
curve obtained from (O and (|18|l . With the parameteri- 
sation for Ns (|18|) we obtain that the critical density, r]c, 
occurs for y/s ~ 10^ GeV. 
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FIG. 1: Particle density as a function of ^/s. Data points are 
from iT^ . The curve was obtained with (I2I18II with y/st = 
10 GeV and a = 1.37. 
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FIG. 2: The inelasticity parameter is' as a function of -^s, 
with K = 1 — XL at relatively low energies and K = 1 — xp 
at energies above the percolation threshold. 

In Fig. 121 we show the ^/s dependence of the inelasticity 
K (equation p7|l ). The behaviour of 1 — xl and 1 — xp 
is also shown in the Figure. From the combination of the 
two curves, K has a maximum at y/s ~ 10^ GeV. 

The development of showers, hadronic and electro- 
magnetic, in cosmic ray physics, is critically dependent 
on the energy carried by the fast particles produced in 
the first hadronic collision, namely the leading particle. 
The shower evolves roughly exponentially, and reaches its 
maximum for a mean depth X (after the first collision). 

Assuming now a simplified branching model, the rela- 
tive shower maximum X can be expressed as: 



X = Xologw [(1 - K)E/Eo] , 



(19) 



where E is the laboratory energy {E ~ ^s), K is the 
inelasticity (as defined in the present percolation model). 
Xq = 60 g/cm^ and Eq = lO'^ eV are effective parameters 
related to the radiation length and to a low energy thresh- 
old for the shower branching, respectively. We have com- 
pared ifT^ to simulations [23 using hadronic interaction 
generators (Sibyll '2l'| and QGSjet based respec- 

tively on the Dual Parton Model and Quark Gluon st ring 
Model), incorporated in the CORSIKA program |2.ll |. 
The comparison is shown in Fig. 13 the parameterisa- 
tion H19|) is, at least for fixed energy, acceptable. 
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FIG. 3: The relative depth of the shower maximum as a 
function of the elasticity parameter 1 — K is shown for Sibyll 
and QGSjet simulations (taken from |23|) and compared to 
the parameterisation of equation Ill9t . 
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FIG. 4: The relative depth of the shower maximum as a func- 
tion of the primary energy. The figure was adapted from , 
superimposing the result of the present percolation model (full 
line). Points are data and dashed lines show predictions of 
QGSJET and SIBYLL for protons and iron. 

The dependence of X^ax = X + Xq on the primary 
energy E is shown in Fig.0| Near the percolation thresh- 
old {E ^ 10^ GeV) there is a clear increase in the slope 
of X{E). In the present percolation model, the change 
in the slope of the X{E) curve can be qualitatively ex- 
plained, in a natural way, through a change in the be- 
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haviour of the inelasticity, due to the effect of percolat- 
ing sea strings above a certain energy. In contrast, this 
change in slope is usually explained by changing the frac- 
tion of heavy nuclei in cosmic rays (see for instance 131): 
this fraction would be higher below the "kink" region (the 
region of fast growth of K), while above it the fraction 
of protons would rise (the region of decrease of K). In 
percolation models, a natural explanation arises without 
requiring a composition change. In fact, Fig. 01 shows 
that the percolation model line naturally goes from the 
model predictions for iron to those for protons as the en- 
ergy increases. A quantitative description of the data is 
however beyond the scope of the present work, as it would 
imply a dedicated Monte Carlo simulation of proton-air 
interactions including percolation effects. 

Furthermore, at very high energies a discrepancy be- 
tween ground array experiments and fluorescence detec- 
tors is usually quoted 2^. We may note that the exis- 
tence of percolation will basically only affect the energy 
measurement of ground arrays, which relies heavily on 
the Monte Carlo simulation of the shower development. 
In fact, including percolation K is lower and fast secon- 
daries are expected. X^ax will thus increase. This will 
lead, at the atmospheric depth of AGASA (taking into 
account that the acceptance is maximal for relatively in- 
clined showers) both to a larger total number of particles 
and to a larger particle density in the region 600 m away 
from the shower core. This effect could partially explain 
the apparent contradiction between ground and fluores- 
cence experiments at GZK energies. 

Additional work covering in detail the general case of 
AB collisions and the treatment of the high energy be- 
haviour of cross-sections is in preparation 27] . 
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